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SEASOAR and CTD Observations During a COARE Surveys Cruise, 
W9211C, 22 January to 22 February 1993. 


Introduction 

An international Coupled Ocean-Atmosphere Response Experiment 
(COARE) was conducted in the warm-pool region of the western equatorial 
Pacific Ocean over a four-month period from November 1992 through February 
1993 (Webster and Lukas, 1992). Most of the oceanographic and meteorological 
observations were concentrated in the Intensive Flux Array (IFA) centered at 1° 
45'S, 156°00'E. As part of this experiment, the R/V Wecoma conducted three 
survey cruises on the R/V Wecoma; each cruise included measurements of the 
temperature, salinity and velocity distribution in the upper 300 m of the ocean, 
and continuous meteorological measurements of wind, air temperature, 
humidity, etc. Most of these measurements were along a butterfly pattern that 
was sampled repeatedly during the three COARE Surveys cruises, W9211A and 
W9211B, and W9211C. 


Coordinates of the Standard Butterfly Pattern were chosen to measure 
zonal and meridional gradients across the center of the IFA, spanning the 
profiling current meter array, while avoiding moorings and stationary ships 
without frequent deviations from our track (Figure 1). The standard coordinates 


of the butterfly apexes are: 
SBN 1°14'S 156°06'E 
SBS 2°26'S 156°06'E 
SBW 1°50'S 155°30'E 
SBE 1°50'S 156°42'E, 


and sampling was done sequentially along the track joining these four points, i.e. 
along a meridional section (N2S) from SBN to SBS, a diagonal section (S2W) 
from SBS to SBW, a zonal section (W2E) from SBW to SBE, and a diagonal 
section (E2N) from SBE to SBN to complete the pattern. Along this track, we 
measured the upper ocean temperature and salinity by means of a towed 
undulating Seasoar vehicle (Figure 2) equipped with a SeaBird CTD system, 
while underway at 7-8 knots. CTD casts were made at the beginning and end of 
each tow, primarily to check calibration of the Seasoar sensors. Water velocity 
along the ship's track was measured by means of the ship-borne acoustic 
Doppler current profiler. 


This report summarizes the Seasoar and CTD observations from 
Wecoma's third COARE Surveys cruise, W9211C. It also provides a cruise 
narrative, and a brief description of the data processing procedures. 
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Figure 1. The Standard Butterfly Pattern in relation to the moorings of the 
COARE Intensive Flux Array. 


Figure 2. Sketch of the Seasoar vehicle used during W9211C. 


Table 1. Summary of CTD stations during W9211C. 


Date Time Station Latitude Longitude Wind Wind Atmos. 


(UTC) No. ir.(T) Spd. (kts) P. (mbar) 
26Jan 1439 1 03°07.4'N 156°07.1'E 260 10 1005.8 
27Jan 0725 2 01°14.0'S 156°06.0' 310 8 1006.8 
29Jan 0540 3 02°27.2' 156°07.8' 300 9 1007.2 
29Jan 0750 4 02°22.0'  156°02.0 320 11 1008.3 
3Feb 2035 5 01°513'  155°30.3' 340 20 1010.8 
3Feb 2159 6 1°51.2'  155°34.1' 350 19 1012.3 
7Feb 0003 7 ~ 01°172'  156°09.6' 340 10 1010.7 
7Feb 1930 8  01°14.0'  156°06.0' 340 12 1008.4 
12Feb 2059 9 02°004' 155°40.4' 335 13 1009.6 
12Feb 2320 10 02°003'  155°40.7 330 15 1010.0 
13Feb 0440 11 01°593'  155°43.4' 310 10 1007.2 
16Feb 2130 12 00°087S 156°06.1'E 070 10 1010.3 


Table 2. Instruments and sensors used for CTD, Seasoar, and underway salinity 
sampling, W9211C, and date of most recent manufacturer's pre-cruise calibration. 


System (Instrument) Sensor Pre-Cruise Calibration Date 


CTD (SBE 9/11 plus SN 0256) P 50130 5 Mar 92 


Tl 1364 6 Oct 92 (modified 2 Dec 92) 
T2 1366 6 Oct 92 (modified 2 Dec 92) 
Cl 1018 16 Sept 92 
C2 1021 16 Sept 92 
Seasoar (SBE 9/11 plus SN 2843) 
(Tows 1-6) P 39107 5 Mar 92 
Cl 1030 17 Apr 92 
C2 1041 24 Apr 92 
(Tows 1-3) Tl 1367 27 Mar 92 (modified 2 Dec 92) 
T2 1369 27 Mar 92 (modified 2 Dec 92) 
(Tows 4-6) Tl 997 6 Oct 92 (modified 2 Dec 92) 
T2 1384 6 Oct 92 (modified 2 Dec 92) 
5-m Intake (MIDAS) T 854 10 Aug 90 
Cc 1070 23Sept 92 


T-C sensor pair Pump 
with duct (inside nose) 


Inlet 
out 
SeaSoar body outline 


Figure 3. Schematic of the plumbing of the ducted T/C sensors inside the Seasoar 
vehicle. Primary sensor inlet and outlet ports were on the starboard side of the Seasoar 
nose; secondary sensor ports were on the port side. During Tow 2, both ducts were 
disconnected at the outlet from the pump; during Tow 4, the secondary sensor duct 
was disconnected at the outlet port. 


Top of upper side rail 


sin u = (U - 1.75) /L 
sin d = (D + 1.75) /L 


Figure 4. Schematics of Seasoar wing angle settings. During Tows 1 and 2, the value of 
D was 23/3" and U was 75/g", yielding an up-angle of 23° and a down-angle of 16°. 
During Tows 3 through 6, D was 25/g" and U was 73/g", yielding an up-angle of 22° and 
a down-angle of 17°. 
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Cruise Narrative, W9211C 

Wecoma departed from Guam about 0000 UTC, 22 Jan 1993, and began a 
transit to 3°N, 156°E where we intended to begin a cross-equatorial Seasoar 
section. After arrival on 26 January, we made a pre-tow CTD cast (Station 1, 
Table 1) with an SBE 9/11 plus CTD with dual ducted temperature and 
conductivity sensors (Table 2). Seasoar was deployed for Tow 1 at about 0600 
UTC, 26 Jan at 3°08'N, 156°07'E; the Seasoar vehicle (Figure 2) was equipped 
with an SBE 9/11plus CTD (SN 0256) with dual ducted temperature and 
conductivity sensors (Table 2; Figure 3). Tow 1 began normally, but in about 10 
minutes the vehicle stopped responding to the control signal, and Seasoar was 
recovered on deck about 0845 UTC, 26 Jan. There was no evidence of a leak in 
the hydraulic unit, but it had definitely lost its power (after working successfully 
for more than 800 hours of consecutive towing in W9211A and W9211B). We 
then proceeded directly to the COARE Intensive Flux Array to begin sampling 
along the Standard Butterfly Pattern (Figure 1). 


A new hydraulic unit (SN 011) was installed in the Seasoar vehicle, and 
the pushrods were adjusted to set the wing-angle so that the maximum up-angle 
was 23° and the maximum down-angle was 16° (Figure 4). 


When we arrived at SBN, the northern apex of the Standard Butterfly 
Pattern, we made a pre-tow CTD cast to 500 m (Station 2, Table 1) and then 
deployed Seasoar at about 0845 UTC, 27 January. Tow 2 began southward 
toward SBS along the N2S line and continued along the Standard Butterfly in the 
usual direction (Table 3); for a 12-km portion of the line between SBE and SBN, 
the Seasoar was kept at depths below 30 m to avoid contamination from 
emptying the ship's tanks. The Seasoar vehicle was very difficult to control in the 
presence of strong vertical shear, so the resulting trajectories were quite irregular 
(e.g., Figure 5); we tentatively concluded that the 16° down-angle was too 
shallow. Preliminary temperature and salinity profiles and T-S diagrams 
showed there was a flow-rate problem through both T-C ducts from about the 
beginning of Tow 2, particularly during descent. Seasoar was recovered at SBS at 
about 0500 UTC, 29 January and we found that the outflow from the pumps was 
not connected to the tubes leading to the outlet ports through the nose of the 
vehicle (Figure 3); instead, the outflow from both sensor ducts was about 10 
inches above the intake and into the interior of the Seasoar vehicle during Tow 2. 
CTD Station 3 (Table 1) was completed immediately after recovery. 


After Tow 2, we changed the Seasoar wing angle setting to yield a 
maximum up-angle of 22° and a maximum down-angle of 17° (Figure 4) and 
reinstalled the CTD, checking the continuity of both sensor ducts with running 
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Table 3. Times (UTC) of standard waypoints during Tow 2 of W9211C. 
Positions of these waypoints are: SBN (1°14'S, 156°06'E); SBS (2°26'S, 156°06'E); 
SBW (1°50'S, 155°30'E); SBE (1°50'S, 156°42'E). 


27 January 
28 January 
28 January 
29 January 
29 January 


Begin/End 


0846 


0503 


SBN 
0917 


1930 


SBW 


0143 


SBE 


1207 


Table 4. Times (UTC) of standard waypoints during Tow 3 of W9211C. 


29 January 
30 January 
30 January 
31 January 
31 January 
1 February 
2 February 
2 February 


Begin/End 


0920 


2004 


SBN 


SBS 


SBW 


1646 


0238 


1124 
1925 


SBE 


0225 


1112 


2015 


Table 5. Times (UTC) of standard waypoints during Tow 4 of W9211C. 


3 February 
4 February 
5 February 
5 February 
6 February 
6 February 


Begin/End 


1306 


2315 


SBN 


0604 


1425 


SBS 


1524 


0015 


SBW 


SBE 
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Table 6. Times (UTC) of standard waypoints during Tow 5 of W9211C. 


Begin/End SBN SBS SBW SBE 
7 February 2058 2200 
8 February 0820 1524 
9 February 0114 
9 February 0803 1824, 1921 
10 February 0200 1128 
10 February 1819 
11 February 0452 1207 2135 
12 February 0433 1456 
12 February 2036 


Table 7. Times (UTC) of standard waypoints during Tow 6 of W9211C. 


Begin /End SBN SBS SBW SBE 
13 February 0551 0845 1856 
14 February 0151 1139 1833 
15 February 0456 
15 February 1137 
15 February 2047 


Table 8. Summary of Seasoar Tows, W9211C. Parameter values in the last column (for 
the T/C alignment offset, and for the amplitude and time constant of the thermal 
mass of the conductivity cell) were used for the preliminary (at-sea) data processing. 


Tow Start Time Stop Time Duration Parameters T/C Pair (and Parameter 
No. (hrs) Measured used for 
At-Sea Analysis 

1 01/26/0605 01/26/0840 0 P, T1, C1, T2, C2 T2, C2 (3.25, 0.045, 8.0) 

2 01/27/0846 01/29/0500 50 P, T1*, C1*, T2*,C2* 12, C2 (3.25, 0.045, 8.0) 

3 01/29/0920 02/02/2004 106 P, T1, C1, T2, C2 T2, C2 (3.25, 0.045, 8.0) 

4 02/08/1306 02/06/2315 80 P, T1,C1, T2*,C2* 11, C1 (2.00, 0.045, 5.0) 

5 02/07/2059 02/12/2063 118 P, T1, C1, T2, C2 T2, C2 (2.00, 0.045, 5.0) 

6 02/13/0551 02/15/2045 62 P, T1, C1, T2, C2 T2, C2 (2.00, 0.045, 9.0) 


*Use data from ascending profiles only; outlet of T-C duct was into the interior of the Seasoar 


vehicle. 
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water. Tow 3 was deployed at about 0915 UTC, 29 January at 2°22'S, 156°02'E 
near SBS immediately after CTD Station 4. Tow 3 began northwestward toward 
SBW along the S2W line, and continued along the Standard Butterfly Pattern in 
the usual direction (Table 4). The Seasoar trajectory (Figure 5) was greatly 
improved over the previous tow, presumably because of the 1° change in the 
wing-angle setting. Preliminary T-S diagrams showed that data quality from 
both sensor pairs was also greatly improved. Preliminary time series of the 
temperature difference between T1 and T2 showed both drift and offsets larger 
than expected from SBE specifications (Figure 6). Tow 3 continued for more than 
4 days (Table 4). Seasoar was recovered near SBW at about 2010 UTC, 2 
February and CTD Station 5 (Table 1) was completed immediately after 
recovery. 


After Tow 3, we replaced both primary and secondary temperature 
sensors with sensors without a reed-switch assembly (SN 997 and 1384) to 
reduce the likelihood of damage from vibration, and reterminated the tow cable. 
Wecoma ran at full speed for a few hours to clear the exhaust system. 


Seasoar was deployed near SBW at about 1300 UTC, 3 February 
immediately after CTD Station 6, and Tow 4 proceeded along the Standard 
Butterfly pattern in the usual direction (Table 5). Preliminary salinity data from 
the secondary sensor pair were very noisy during ascending profiles, and we 
suspected a problem with the sensor duct. Preliminary time series of the 
difference between the two temperature sensors showed no evidence of drift or 
offset (Figure 6). The temperature and conductivity data became intermittent 
during the E2N section of 6 February, so the vehicle was recovered near SBN at 
about 2315, 6 Feb and CTD Station 7 was made soon afterward. We found that 
the outlet tube for the secondary sensors had become detached from the nose, so 
the actual outlet from these sensors was inside the vehicle. 


After Tow 4, we reterminated the Seasoar cable again, using a modified 
wire grip instead of the standard tapered plastic cow-tail to reduce cable flexing. 
We also repaired the secondary outlet duct, so that this outlet was once again 
directly through the nose of the Seasoar vehicle. Seasoar was deployed near SBN 
at about 2100, 7 February, immediately after CTD Station 8 (Table 1). Tow 5 
began southward toward SBS and continued along the Standard Butterfly in the 
usual direction for almost five days (Table 6), when it was terminated at about 2° 
on the S2W line to allow the ship to make a high-speed run to clear the exhaust 
system; CTD Station 9 was made immediately after recovery. The Seasoar 
vehicle, instrument and termination were apparently in excellent condition. 
During CTD Station 10, as we were preparing to deploy Seasoar for Tow 6, we 
found a broken strand in the Seasoar cable adjacent to the loop attached for 


crane handling. We therefore cut the cable above this point and reterminated the 
cable. 


Seasoar was deployed for Tow 6 at about 0545 UTC, 13 February at 1°59'S, 
155°43'E on the S2W line, immediately after CTD Station 11. Tow 6 began 
northwestward toward SBW and continued along the Standard Butterfly in the 
usual direction (Table 7) until 1137 UTC, 15 February, when we turned 
northward toward the equator. Seasoar was recovered just south of the equator 
at about 2045 UTC, 15 February, and CTD Station 12 was completed 
immediately afterward. We then ceased operations in the COARE IFA, and 
Wecoma proceeded to Pohnpei. 


In all, we completed six Seasoar tows during W9211C (Table 8), for a total 
towing time of 416 hours. The overall Seasoar sampling included 12 occupations 
each of the N2S and W2E lines (Table 9) and the S2W and E2N lines (Table 10). 


Wecoma arrived in Pohnpei at about 2300 UTC, 17 Feb, to disembark 
some personnel and departed there at about 0600 UTC, 18 Feb for the transit to 
Guam. En route to Guam, we conducted an 8-hour test of Seasoar flight 
characteristics with an optical plankton counter attached to the vehicle. Wecoma 
arrived in Guam at 0400 UTC, 22 Feb. 


Underway measurements were made continuously through most of the 
cruise. These include: Acoustic Doppler Current Profile measurements of water 
velocity relative to the ship and accompanying GPS position data (contact Eric 
Firing et al., Univ. of Hawaii); temperature and salinity of water at 5 m depth 
and 2 m depth (contact Clayton Paulson, Oregon State Univ.); near-surface 
salinity of water pumped from a buoyant hose (contact Gary Lagerloef, SAIC); 
microwave radiometer estimates of surface salinity (contact Gary Lagerloef, 
SAIC); and a broad spectrum of meteorological observations (contact Clayton 
Paulson, OSU) including sonic inertial dissipation (contact Jim Edson, WHO). 


Members of the scientific party included Marc Willis, Tim Holt (both 
Wecoma Marine Technicians), Adriana Huyer, Robert L. Smith, Fred Bahr, Pip 
Courbois, Jane Fleischbein (all from Oregon State University), Eric Firing, Fred 
Bingham, Tung Le, Dail Rowe, and Joanna Muench, (all from University of 
Hawaii), Clay Wilson (SAIC) and Jonas Aleksa (University of Massachusetts). 
Additional persons on the Pohnpei-Guam transit were Mike Hill (Oregon State 
University), Meng Zhou, and Walter Nordhausen (both from Scripps Institution 
of Oceanography). 
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Table 9. Times (UTC) of meridional and zonal sections of the Standard Butterfly 
pattern. All N2S sections were southward along 156°06'E from SBN 

(1°14'S) to SBS (2°26'S), and all W2E sections were eastward along 1°50'S 

from SBW (155°30'E) to SBE (156°42'E). Number in parentheses indicates intake of 
preferred sensor T/C sensor pair was on port (1) or starboard (0) side of Seasoar) 


N2S (SBN to SBS) along 156°06'E 


0917, 27 Jan to 1827, 27 Jan* (0) 
1930, 28 Jan to 0430, 29 Jan*(0) 
0945, 30 Jan to 1926, 30 Jan (0) 
1827, 31 Jan to 0358, 1 Feb (0) 
0320, 2 Feb to 1221, 2 Feb (0) 
0604, 4 Feb to 1524, 4 Feb (0) 
1425, 5 Feb to 0015, 6 Feb (0) 
2200, 7 Feb to 0820, 8 Feb (1) 
0803 to 1824, 9 Feb! (1) 
1819, 10 Feb to 0452, 11 Feb (1) 
0433 to 1456, 12 Feb (1) 
0151 to 1139, 14 Feb (1) 
* Data from ascending profiles only. 


W2E (SBW to SBE) along 1°50'S 


0143, 28 Jan to 1207, 28 Jan* (0) 
1646, 29 Jan to 0225, 30 Jan** (0) 
- 0238 to 1112, 31 Jan (0) 
1124 to 2015, 1 Feb (0) 
1410 to 2304, 3 Feb (0) 
2230, 4 Feb to 0719, 5 Feb (0) 
0747 to 1641, 6 Feb (0) 
1524, 8 Feb to 0114, 9 Feb (1) 
0200 to 1128, 10 Feb (1) 
1207 to 2135, 11 Feb (1) 
0845 to 1856, 13 Feb (1) 
1833, 14 Feb to 0456, 15 Feb (1) 


** Includes a 6 second data gap from 20:52:25 through 20:53:30 UTC, 29 January. 


1 Section extends past SBS to 2°29'S. 


Table 10. Times (UTC) of diagonal sections of the Standard Butterfly pattern: 


S2W between SBS (2°26'S, 156°06'E) and SBW (1°50'S, 155°30"E); and E2N 
between SBE (1°50'S, 156°42'E) and SBN (1°14'S, 156°06'E). During most E2N 
sections Seasoar was kept at maximum depth for about 12 km. 


S2W (SBS to SBW) 


1827, 27 Jan to 0142, 28 Jan* (0) 
1035, 29 Jan to 1646, 29 Jan (0) 
1926, 30 Jan to 0238, 31 Jan (0) 
0358, 1 Feb to 1124, 1 Feb (0) 
1221, 2 Feb to 1925, 2 Feb (0) 
1524 to 2230, 4 Feb (partial) (0) 
0015 to 0747, 6 Feb (0) 
0820 to 1524, 8 Feb (1) 
1824, 9 Feb to 0200, 10 Feb (1) 
0452 to 1207, 11 Feb (1) 


1456 to 2036, 12 Feb (partial) (1) 
0641 to 0845, 13 Feb (partial) (1) 


1139 to 1833, 14 Feb (1) 
*Data from ascending profiles only. 


E2N (SBE to SBN) 


1207, 28 Jan to 1930, 28 Jan* (0) 
0225 to 0945, 30 Jan (0) 
1112 to 1827, 31 Jan (0) 

2015, 1 Feb to 0320, 2 Feb (0) 
2304, 3 Feb to 0604, 4 Feb (0) 
0719 to 1425, 5 Feb (0) 

1641 to 2315, 6 Feb (partial)** (0) 
0114 to 0803, 9 Feb (1) 
1128 to 1819, 10 Feb (1) 
2135, 11 Feb to 0433, 12 Feb (1) 


1856, 13 Feb to 0151, 14 Feb (1) 
0456 to 1137, 15 Feb (1) 


** Includes numberous short gaps after 22:06:27 UTC, 6 February. 
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CTD Data Acquisition, Calibration and Data Processing 


All CTD/rosette casts were made with an SBE 9/11-plus CTD system 
equipped with dual ducted temperature and conductivity sensors (Table 2). CTD 
casts were made primarily to monitor the calibration of the Seasoar data, and 
were therefore made before and after each Seasoar tow, with as little delay as 
possible. Maximum sampling depth was about 500 m. Raw 24 Hz CTD data 
were acquired on an IBM compatible PC using the SEASAVE module of 
SEASOFT version 4.015 (Anon., 1992); temperature and conductivity data were 
recorded from both pumped sensor ducts. At each station a few salinity samples 
were collected for in situ calibration of the conductivity sensors; CTD values at 
the sample depth (calculated from the most recent manufacturer's pre-cruise 
calibration) were recorded both by pressing the F5 key at the time of rosette 
firing and manually on the station log sheets. Samples were analyzed on a 
Guildline Autosal salinometer that was standardized with LAPSO Standard 
Water P-119 at the beginning and end of each batch of about 36 samples. Sample 
salinity values were in essential agreement with the calculated CTD salinity 
values from both sensor pairs (average AS1=0.002, with standard deviation of 
0.004, 36 samples: average AS2=0.002, with standard deviation of 0.005, 33 
samples), and hence no conductivity correction was applied. 


CTD data were processed on an IBM-compatible PC using applicable 

_ SEASOFT modules. Since there was no significant difference between the data 
from the two sensor pairs, we fully processed data from the primary sensors 
only. The DATCNV module of SEASOFT was used with the pre-cruise 
calibration constants to calculate 24 Hz values of pressure, temperature and 
conductivity from the raw frequencies. When necessary, the output data file was 
edited to remove any spikes and any values inadvertently recorded before the 
pressure minimum at the beginning of the cast. The CELLTM module was used 
to correct for the thermal mass of the conductivity cell, assumed to have a 
thermal anomaly amplitude of 0.03 and a time constant of 9 seconds. Ascending 
portions of the 24-Hz data file were removed by LOOPEDIT with the minimum 
velocity set to 0.0 m/s. The remaining data were averaged to 1 dbar values using 
BINAVG. The final processed data files consist of 1 dbar values of pressure, 
temperature and conductivity. These processed data files were transferred to a 
SUN computer where we used standard algorithms (Fofonoff and Millard, 1983) 
to calculate salinity, potential temperature, density anomaly (sigma-theta), 
specific volume anomaly, and geopotential anomaly (dynamic height). 
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Seasoar Data Acquisition and Preliminary Processing 


Raw 24 Hz CTD data from the Seasoar vehicle and GPS position and time data 
were acquired by an IBM compatible PC, which also set flags to indicate missing GPS 
data and to mark the collection of a salinity sample. The raw data were simultaneously 
recorded on optical disk by PC and on a Sun Sparc workstation. The PC displayed 
time series of subsampled temperature (both sensors), conductivity (both sensors) and 
pressure in real time; it also displayed accumulated temperature data for 6-8 hours as 
a vertical section (color raster). One-second averages of ship's position, CTD 
temperature (both sensors), conductivity (both sensors), salinity (both sensor pairs), 
and pressure were calculated on the Sparc workstation, using the most recent 
manufacturer's calibration (Table 2). For each tow, the preliminary salinity for each 
sensor pair was calculated using a fixed offset between temperature and salinity, and a 
fixed value for the amplitude and time constant of the thermal mass of the 
conductivity cell, but these parameters were changed from one tow to another (Table 
8). Time-series and vertical profile plots of the one-second data were made at the end 
of each hour. The 1-second preliminary data were used to calculate the average 
temperature and salinity data in bins extending 3 km in the horizontal and 2 dbar in 
the vertical, and these gridded values were used to plot vertical sections for each leg of 
the Standard Butterfly pattern. 


Seasoar Conductivity Calibration 


Salinity samples were collected once per hour from a throughflow system in 
Wecoma's wetlab from 0000 UTC, 27 January until 2000 UTC, 15 February 1993. This 
system pumps water from the seachest at a depth of 5 m in the ship's hull, through a 
tank containing SBE temperature and conductivity sensors; samples are drawn from a 
point just beyond this tank. The 120 ml glass sample bottles were rinsed three times 
before filling, and closed with screw-on plastic caps with conical polyethylene liners. 
Samples were further sealed by wrapping parafilm around the base of the cap. 
Samples were analyzed at sea on an Autosal salinometer, usually within 2-3 days after 
collection; the salinometer was standardized with LAPSO Standard Water P-119 at the 
beginning and end of each batch of about 24 samples. Time series of these hourly 
salinity samples and time series of the preliminary Seasoar data from the 3-7 m depth 
range (Figure 7) show very similar variations, especially during Tows 4-6. 


For a quantitative comparison between the salinity samples and the Seasoar 
data, we selected Seasoar values that were both within 7 minutes of the time of the 
salinity sample and within a depth range of 3.0 to 5.5 m. For each salinity sample, we 
calculated a bottle conductivity using the measured salinity and the temperature from 
each Seasoar sensor duct, and then compared this sample conductivity to the directly 
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Figure 8 (a). Time series of salinity differences between the 5-m samples and the 
matching corrected Seasoar data, for the primary T/C sensor pair during Tows 2 
and 3, and for both primary and secondary sensors during Tow 4 of W9211C. 
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measured conductivity from the same sensor duct; a very few pairs with very large 
differences (only 5 out of a total of 329) were eliminated from the comparison. 
Assuming, as usual, that the measured conductivity should be corrected by a 
multiplier alone, we calculated the slope (m) of the zero-intercept regression line 
between the measured conductivity and the sample conductivity separately for each 
sensor pair and for each Seasoar tow. 


Between Tows 2 and 3, we found no significant difference in the multiplier for 
the primary sensor pair, and therefore we pooled the data, and used the same 
correction factor (k) for both tows (Table 11). Since the secondary temperature sensor 
malfunctioned during Tows 2 and 3, we used the primary temperature data to 
determine the multiplier for secondary conductivity for these tows. Both primary and 
secondary sensors were replaced after Tow 3. 


Between Tows 4, 5 and 6, we found no significant difference in the multiplier 
(m) for either the primary or secondary conductivity data, and therefore we adopted 
the values determined from the pooled data as the correction factors (Table 11). 


Applying these multipliers to the Seasoar conductivity data before 
recalculating salinity allows us to compare the corrected Seasoar salinity values from 
both primary and secondary sensor ducts to the sample salinity (Figure 8, Table 11). 
The time series of the differences (Figure 8) show reasonable agreement between 
sample and near-surface Seasoar data for the primary sensor pairs during Tows 2 and 
3, and for both sensor pairs during Tows 4, 5 and 6. In general, largest sample-Seasoar 
differences occur when the surface layer is stratified, so the standard deviations of the 
salinity differences (Table 11) primarily reflect sampling errors rather than 
instrumental noise. 


Table 11. Conductivity multipliers (m1 and m2) for primary and secondary sensors, 
determined (separately for each tow) from comparison of near-surface Seasoar data with 5-m 
intake samples, and the conductivity correction factors (k1 and k2) adopted for reprocessing 
the Seasoar conductivity data. Also shown are the average and standard deviations of the 
salinity differences between the sample values and the corrected Seasoar data. 


Average Std. Dev. 


mi m2 kl k2 S1 $2 $1 $2 

= - 1.000433 1.000449 — - - - 

1.000470 —- 1.000433 1.000449 0.001 - 0.007 — 
= 1.000433 1.000449 0.000 — 0.009 — 


1.000388 1.000460 1.000368 1.000418 0.001 0.002 0.004 0.004 
1.000354 1.000396 1.000368 1.000418 -0.001 -0.001 0.007 0.008 
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Figure 8 (b). Time series of salinity differences between the 5-m samples and the 
matching corrected Seasoar data, for both primary and secondary sensor pairs 
during Tows 5 and 6 of W9211C. 


Post-processing of Seasoar Data 


Background Salinity data derived from SeaBird ducted temperature and 
conductivity sensors are subject to errors from three separate sources (Larson, 1992): 
(1) poor alignment of the 24 Hz temperature and conductivity data, (2) poor 
compensation for the transfer of heat between the mantle of the conductivity cell 
and the water flowing through it, and (3) mismatch of the effective time constants of 
the temperature and conductivity measurements. These sources of error are 
minimized in a normal SeaBird CTD, by pumping the water through the ducted 
pair at a fixed rate, by software algorithms for changing the T/C alignment and the 
thermal mass parameters, and by matching the thermistor to the effective time 
constant of the conductivity measurement through the fixed-length duct. With 
sensors mounted inside a towed Seasoar vehicle, these error sources need to 
considered anew. 


While at sea, we noticed that there was often a significant difference in data 
quality between ascending and descending profiles, with salinity data from 
descending profiles generally noisier than salinity data from ascending profiles; 
descending profiles frequently showed a large and variable positive salinity bias 
(e.g. Figure 9). These systematic up/down differences were observed in data from 
both sensor pairs, though sometimes more pronounced in one of the sensor pairs. 
We therefore suspected that there might be variations in the flow-rate of water 
through the pumped T/C sensor ducts, leading to variations in the optimum 
alignment offset between the 24-Hz temperature and conductivity data. 


To check this hypothesis, we applied the general procedures outlined by 
Nordeen Larson (1992) to a segment of raw data from the hour beginning at 0000 
UTC, 9 February. We first reprocessed the data from both sensor pairs without 
applying any T/C offset, retaining 24-Hz output; data from the first ascending 
profile (about 0001 to 0006 UTC) and the first descending profile (0006 to 0010 UTC) 
were read into separate files. These served as input to a program which realigns the 
temperature and conductivity data from both sensor pairs by a specified offset and 
recalculates salinity; for each scan it also determines a "median salinity", S,,, using a 
centered nine-point median filter, and a “spike salinity", S', as the anomaly from the 
median salinity. The time integral of this spike salinity, normalized by the sign of 
the temperature gradient, 


S=2 (S',/ sgn(Tp41- n-1))s 


changes sign at the optimum offset value (Larson, 1992). Calculations were repeated 
for a range of offsets, and spike salinity profiles were plotted for each (e.g., Figure 
10). For these profiles, we found that the optimum alignment offset, £, during 
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ascent was essentially the same as those for ducted sensors in a normal profiling 
SBE9/11 plus CTD, i.e. 0.0 scans for primary sensors and 1.75 scans for secondary 
sensors. (The SBE 11/plus deck unit previously applies a T/C alignment offset of 
1.75 scans to data from primary sensors but not secondary sensors). Optimum 
values of the offset during descent were larger: 1.75 scans for primary sensors, and 
2.75 scans for secondary sensors (Figures 10, 11). 


We independently calculated lagged correlation coefficients between the 
temperature and conductivity data from each sensor pair, using first differencing to 
remove low frequency trends in both, and obtained the same results, and found that 
the lag of maximum correlation was in essential agreement with the alignment 
offset determined by Larson's spike salinity method. Alternating between two 
values of the alignment offset (one for descent, one for ascent) significantly 
improved the processed salinity data, but systematic up/down differences 
remained during portions of some Seasoar tows (e.g. Figure 12). 


Further study of the 24-Hz Seasoar CTD data has shown that the optimum 
T/C alignment offset, &, varies on several different time scales: within each 
undulation, from section to section, gradually over the length of a tow, and from 
one tow to another (e.g., Figure 13). In most cases, the optimum alignment offset 
during ascent is near the default values of 0.0 scans for T1/C1 and 1.75 scans for 
T2/C2, while the optimum offset during descent is longer. We attribute these 
variations in optimum offset to reductions in the rate of pumped flow through the 
sensors, and believe these variations in flow rate are the result of dynamic pressure 
differences, partly between the inside and outside of the vehicle, and partly along 
the exterior of the vehicle nose (though inlet and outlet ports are separated by only a 
few centimeters). Possible sources of such pressure gradients include high rates of 
ascent and descent (sometimes >3 m/s, superimposed on a horizontal tow speed of 
4 m/s), and local small perturbations of the flow field around the vehicle, associated 
perhaps with a persistent roll angle or strong cross-currents. Whatever their source, 
the variations in offset alignment are a major source of potential error in calculating 
salinity, and must be taken into account in reprocessing the data. After some trial 
and error, we chose to determine the optimum lag in three depth ranges (excluding 
the surface mixed layer) for each ascending and descending profile; calculating lags 
for successive 10-second segments proved to be less satisfactory in minimizing 
salinity noise and bias. 


Still further study of the preliminary Seasoar data showed that systematic 
up/down salinity differences were particularly strong near sharp thermoclines, 
even though we had attempted to correct for the thermal mass of the conductivity 
cell (Lueck, 1990), using constant values for the time constant (t) and amplitude (a) 
parameters in the recursive algorithm provided by SeaBird: 
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Figure 10. Time-series of 24-Hz “median salinity” and “spike salinity” from ascending (left) and 
descending (right) segments of Seasoar data between 0001 and 0010 UTC, 9 February, 
calculated for three different values of T/C alignment offsets 1. The “median salinity” is the 
median of 9 scans, and the “spike salinity” is the anomaly from this median. 
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Figure 11. Values of the integral of spike salinity, ISS, for ascending and 
descending segments of data from both T/C sensor pairs, during 0001 
and 0010 UTC, 9 February 1993; note that the value of the optimum offset 
(at ISS = 0) is larger during descent. 
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Figure 12. T-S characteristics for Seasoar data from two different hours, 
processed by two different methods: with constant values of the T/C 
alignment offset (left, §, = 2.0 scans), and with alternating values of the offset 
(right; ) = 2.5 scans during ascent and & = 5.0 scans during descent). In both 
methods, we used constant values of the thermal mass parameters (& = 0.045 
and t = 9 sec). 
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Figure 13. Time series of the lag at maximum cross-correlation between 24-Hz temperature and 
conductivity data (detrended by first-differencing), for both primary and secondary sensor 
pairs during Seasoar Tow 6. Correlations were calculated for three different layers (50-120 dbar, 
plus signs; 120-180 dbar, squares; 180-240 dbar, triangles), and separately for ascending and 
descending profiles; times of turns in the ship's track are indicated by vertical lines. 
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dt = temperature - previous temperature 
ctm = - b* previous ctm + a * dcdt * dt 
corrected conductivity = conductivity + ctm 


where a =2 a / (0.0417 B + 2), dcdt = 0.1 + 0.0006 (temperature - 20), B = 1/t and b= 
1-2a/a. Constant values of « and t had been chosen by two different methods: by 
comparing CTD casts with Seasoar deployment and recovery profiles, and by 
minimizing differences between ascending and descending salinity data near the 
top of the thermocline for some data segments with exceptionally steep gradients at 
the bottom of the surface mixed layer; both methods indicated the best constant 
values were o = 0.04 for the thermal anomaly amplitude and t = 10 sec for the time 
constant. The association of residual salinity bias with sharp thermoclines 
suggested that variations in flow rate through the ducted sensors effectively altered 
the rate of heat transferred from the walls of the conductivity cell. Experimenting 
with different values of indicated that larger values of are appropriate for larger 
values of the alignment offset (Figure 14). 


Finally, we considered the effect of elongating the intake duct (from about 1 
cm on a standard SeaBird T/C ducted pair to about 5 cm) by adding Tygon tubing 
between the sensor duct and inlet in the Seasoar nose. In principle, elongating the 
intake could slow the response time of the conductivity measurement through 
turbulent mixing within the duct (Nordeen Larson, personal communication). We 


_ compared conductivity spectra from conventional CTD casts with those from the 


Seasoar deployment and recovery profiles, and found no significant difference: both 
were similar to the temperature spectra. We conclude that any change in sensor time 
constant due to intake elongation is negligible, at least with the present choice of 
tubing. 


Procedures The first step in reprocessing is to incorporate temperature and 
conductivity corrections determined from in-situ calibration or post-cruise 
calibration in a new configuration file for each tow. Except for temperature sensor 
#1369, which clearly malfunctioned during the cruise, the post-cruise temperature 
calibrations of the Seasoar sensors showed little or no drift from the pre-cruise 
calibrations, so no temperature corrections were applied. The conductivity 
correction factors determined from the in situ calibration data (Table 11) were 
incorporated in the configuration files for each tow. 


The next step is to compute lagged correlations between temperature and 
conductivity for each sensor pair, separately for ascending and descending profiles, 
and separately for three depth ranges: 50 to 120 dbar, 120 to 180 dbar, and 180 to 240 
dbar, provided the segment contains at least 72 scans. Cross-correlations are 
calculated after detrending both temperature and conductivity by first-differencing 
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the 24-Hz data. Correlations are calculated for +12 lags; the maximum correlation is 
almost always 20.85. The fractional value of the lag at maximum correlation is 
determined by fitting a parabola to the cross-correlation values. The resulting time 
series of the optimum primary and secondary alignment offsets (€, and & ) for each 
tow are shown in Appendix A. Lag values greater than 12 and other outliers 
(obtained occasionally for data segments lasting <10 seconds) were not used in 
processing the data. Very large lags were obtained for many descending profiles of 
Tow 2) when the pump outlets were not properly connected, and hence data from 
descending profiles from this tow were discarded. 


The edited values of the alignment offset were applied sequentially in 
reprocessing the 24-Hz T/C data. To reprocess data from depths shallower than 50 
m, we used the € value determined from the preceding 120 to 50 dbar layer; for data 
deeper than 240 m, we used the € value determined from the preceding 180 to 240 
dbar layer. Short segments with unreasonably large lags were processed with the 
lag obtained for the succeeding data segment. 


To correct the 24 Hz conductivity data for the thermal mass of the 
conductivity cell, we used the standard algorithm with a fixed value for the thermal 
anomaly time constant (t = 10 sec), and variable values for the thermal anomaly 
amplitude depending on the alignment offset: 


Q,= 0.03 if Ei <0 
Qa = 0.03 + 0.03(§; /2.75) if & >0 
ay = 0.03 if E<1.75 


Gy = 0.03 + 0.03(E,-1.75)/2.75)) if E, >0, 


where the subscripts denote primary or secondary sensors. The corrected and 
realigned 24 Hz temperature and conductivity data are used to calculate 24-Hz 
salinity, and these are averaged to yield 1-second averages stored in hourly files. 
These procedures greatly reduce the systematic salinity differences (both bias and 
noise) between ascending and descending profiles (e.g. Figure 14). However, 
salinity data from descending profiles remains somewhat noisy, while the 
ascending data is almost always free of significant noise. 


The reprocessed data from both sensor pairs were plotted to determine which 
pair provided the better data for each Seasoar tow. During Tow 2, with both sensor 
ducts disconnected from the outlet, there is evidence of stalling in the pumped 
sensor ducts during descent, but data from ascending profiles is apparently of good 
quality. For both Tows 2 and 3, with the malfunctioning secondary temperature 
sensor, data from the primary sensors is clearly preferable. During Tow 4, the 


secondary sensor duct was disconnected from the outlet, and data from the primary 
sensors is clearly preferable. For Tows 5 and 6, the data from the two sensor pairs 
was of essentially the same quality, with only very subtle reasons to prefer the 
secondary sensors over the primary sensors. 


Data Presentation 


Successive hourly files of the reprocessed one-second average data were 
joined and clipped to yield a single data file for each section of the Standard 
Butterfly Pattern (Tables 9 and 10). Final processed data files contain unfiltered GPS 
latitude and longitude; pressure; temperature, salinity and sigma-t from the better 
sensor pair; date and time; and an integer representing flags (to indicate collection 
of a water sample from 5-m intake (thousands digit set to 1), missing GPS data filled 
by linear interpolation (tens digit set to 1), and to indicate port or starboard intake 
for the T/C sensor pair (ones digit set to 1 or 0, respectively)). We present 
consecutive figures of the Seasoar trajectory (time series of pressure, latitude and 
longitude) along each section. We also present summary figures of all of the 1- 
second data for each of the four standard sections as follows: ensembles of 
temperature profiles (both ascending and descending, except for Tow 2), salinity 
profiles (ascending profiles only), and T-S diagrams (for ascending profiles only). 
Vertical distributions of the temperature, salinity and sigma-t along each section 
were plotted using Don Denbo's PPlus program with a vertical grid spacing of 2 
dbar and a horizontal spacing of 1 nm, and with a value of CAY= 5 for the 
smoothing parameter (combined spline and laplacian filter). For the temperature 
sections, we used both ascending and descending data, for all tows except Tow 2 
where we used ascending data only. For the salinity and sigma-t sections, we used 
only ascending data for all tows. 


CTD/Seasoar Comparison 


T-S diagrams for the beginning and end of each Seasoar Tow are shown in 
Appendix B. Each diagram shows the T-S curve from both the conventional CTD 
cast and the preferred Seasoar sensors during Seasoar deployment or recovery. 
Seasoar deployment profiles are generally noisier than either the CTD profiles or 
Seasoar recovery profiles, probably because the Seasoar vehicle is tilted nose- 
upward during both deployment and recovery; since the ship is moving very 
slowly, observations during deployment are sometimes within the turbulent wake 
of the descending vehicle. 
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CTD DATA 


For each station, we present a plot of the vertical temperature, salinity, and 
sigma-t profiles, and a listing of the observed and derived variables at standard 
pressures. Header data includes the CTD Station Number, Latitude (in degrees 
and minutes; a negative degree integer indicates southern hemisphere), 
Longitude (degrees and minutes East), Date and Time (UTC), and Bottom Depth 
(in meters). 
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